Charlottesville Region
cville_tracts <- readRDS("Cville_Tree/cville_tracts.RDS")
cville_blkgps <- readRDS("Cville_Tree/cville_blkgps.RDS")
cville_blocks <- readRDS("Cville_Tree/cville_blocks.RDS")
Tracts
cvl_tracts <- read.csv("Cville_Tree/nlcd_tree_cville_tracts.csv")%>%
dplyr::select(-c(5,7,8))
cvl_tracts <- cvl_tracts[order(cvl_tracts$TRACTCE), ]
kable(cvl_tracts)%>%
kable_styling(font_size = 10)
|
|
STATEFP
|
COUNTYFP
|
TRACTCE
|
GEOID
|
NAMELSAD
|
tree_can
|
|
19
|
51
|
540
|
201
|
51540000201
|
Census Tract 2.01
|
32.76820
|
|
41
|
51
|
540
|
202
|
51540000202
|
Census Tract 2.02
|
11.58013
|
|
42
|
51
|
540
|
302
|
51540000302
|
Census Tract 3.02
|
25.91021
|
|
43
|
51
|
540
|
401
|
51540000401
|
Census Tract 4.01
|
34.10494
|
|
45
|
51
|
540
|
402
|
51540000402
|
Census Tract 4.02
|
22.95314
|
|
44
|
51
|
540
|
501
|
51540000501
|
Census Tract 5.01
|
34.68436
|
|
21
|
51
|
540
|
502
|
51540000502
|
Census Tract 5.02
|
59.82563
|
|
46
|
51
|
540
|
600
|
51540000600
|
Census Tract 6
|
33.04271
|
|
47
|
51
|
540
|
700
|
51540000700
|
Census Tract 7
|
53.04950
|
|
48
|
51
|
540
|
800
|
51540000800
|
Census Tract 8
|
47.61878
|
|
49
|
51
|
540
|
900
|
51540000900
|
Census Tract 9
|
39.40146
|
|
14
|
51
|
540
|
1000
|
51540001000
|
Census Tract 10
|
24.29448
|
|
37
|
51
|
3
|
10100
|
51003010100
|
Census Tract 101
|
71.90352
|
|
9
|
51
|
3
|
10201
|
51003010201
|
Census Tract 102.01
|
59.63794
|
|
8
|
51
|
3
|
10202
|
51003010202
|
Census Tract 102.02
|
54.97607
|
|
32
|
51
|
3
|
10300
|
51003010300
|
Census Tract 103
|
57.39112
|
|
10
|
51
|
3
|
10401
|
51003010401
|
Census Tract 104.01
|
61.46941
|
|
11
|
51
|
3
|
10402
|
51003010402
|
Census Tract 104.02
|
61.50637
|
|
34
|
51
|
3
|
10500
|
51003010500
|
Census Tract 105
|
63.73709
|
|
50
|
51
|
3
|
10601
|
51003010601
|
Census Tract 106.01
|
44.70170
|
|
30
|
51
|
3
|
10602
|
51003010602
|
Census Tract 106.02
|
50.20944
|
|
33
|
51
|
3
|
10700
|
51003010700
|
Census Tract 107
|
29.62388
|
|
39
|
51
|
3
|
10800
|
51003010800
|
Census Tract 108
|
59.78755
|
|
40
|
51
|
3
|
10901
|
51003010901
|
Census Tract 109.01
|
33.77119
|
|
24
|
51
|
3
|
10902
|
51003010902
|
Census Tract 109.02
|
71.46157
|
|
31
|
51
|
3
|
10903
|
51003010903
|
Census Tract 109.03
|
17.84085
|
|
35
|
51
|
3
|
11000
|
51003011000
|
Census Tract 110
|
66.04177
|
|
38
|
51
|
3
|
11100
|
51003011100
|
Census Tract 111
|
47.54826
|
|
6
|
51
|
3
|
11201
|
51003011201
|
Census Tract 112.01
|
72.63061
|
|
7
|
51
|
3
|
11202
|
51003011202
|
Census Tract 112.02
|
68.84397
|
|
12
|
51
|
3
|
11301
|
51003011301
|
Census Tract 113.01
|
67.37963
|
|
13
|
51
|
3
|
11302
|
51003011302
|
Census Tract 113.02
|
61.29524
|
|
5
|
51
|
3
|
11303
|
51003011303
|
Census Tract 113.03
|
48.20600
|
|
36
|
51
|
3
|
11400
|
51003011400
|
Census Tract 114
|
67.93056
|
|
29
|
51
|
65
|
20101
|
51065020101
|
Census Tract 201.01
|
61.97026
|
|
28
|
51
|
65
|
20102
|
51065020102
|
Census Tract 201.02
|
64.89842
|
|
23
|
51
|
65
|
20200
|
51065020200
|
Census Tract 202
|
61.53652
|
|
22
|
51
|
65
|
20300
|
51065020300
|
Census Tract 203
|
65.79589
|
|
20
|
51
|
79
|
30101
|
51079030101
|
Census Tract 301.01
|
67.35107
|
|
4
|
51
|
79
|
30102
|
51079030102
|
Census Tract 301.02
|
69.19865
|
|
3
|
51
|
79
|
30200
|
51079030200
|
Census Tract 302
|
47.86960
|
|
17
|
51
|
125
|
950100
|
51125950100
|
Census Tract 9501
|
76.19516
|
|
25
|
51
|
109
|
950100
|
51109950100
|
Census Tract 9501
|
56.51193
|
|
16
|
51
|
125
|
950200
|
51125950200
|
Census Tract 9502
|
78.27371
|
|
15
|
51
|
109
|
950201
|
51109950201
|
Census Tract 9502.01
|
55.55105
|
|
1
|
51
|
109
|
950202
|
51109950202
|
Census Tract 9502.02
|
59.46266
|
|
2
|
51
|
109
|
950300
|
51109950300
|
Census Tract 9503
|
56.50907
|
|
18
|
51
|
125
|
950300
|
51125950300
|
Census Tract 9503
|
74.95332
|
|
27
|
51
|
109
|
950400
|
51109950400
|
Census Tract 9504
|
60.34939
|
|
26
|
51
|
109
|
950500
|
51109950500
|
Census Tract 9505
|
61.99446
|
cville_tracts$tree_can <- cvl_tracts$tree_can
pal <- colorNumeric("plasma", reverse = TRUE, domain = cvl_tracts$tree_can)
cvl_tracts$tree_can <- as.numeric(cvl_tracts$tree_can)
cville_tracts$COUNTYFP <- as.numeric(cville_tracts$COUNTYFP)
m <- leaflet()%>%
addTiles()%>%
addPolygons(data = cville_tracts,
fillColor = ~pal(tree_can),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(weight = 2, fillOpacity = 0.8, bringToFront = T),
popup = paste0("FIPS Code: ", cvl_tracts$GEOID, "<br>",
"Tree Canopy: ", cvl_tracts$tree_can,"<br>",
"Tract: ", cvl_tracts$NAMELSAD)) %>%
addLegend("bottomright", pal = pal, values = cvl_tracts$tree_can,
title = "Tree Canopy", opacity = 0.7)
m
Blockgroups
cvl_blkgps <- read.csv("Cville_Tree/nlcd_tree_cville_blkgps.csv")%>%
as.data.frame()
cvl_blkgps <- cvl_blkgps[order(cvl_blkgps$TRACTCE), ]
kable(head(cvl_blkgps))
|
|
STATEFP
|
COUNTYFP
|
TRACTCE
|
BLKGRPCE
|
GEOID
|
NAMELSAD
|
tree_can
|
|
54
|
51
|
540
|
201
|
3
|
5.154e+11
|
Block Group 3
|
39.102459
|
|
55
|
51
|
540
|
201
|
1
|
5.154e+11
|
Block Group 1
|
43.246032
|
|
89
|
51
|
540
|
201
|
2
|
5.154e+11
|
Block Group 2
|
15.508475
|
|
6
|
51
|
540
|
202
|
2
|
5.154e+11
|
Block Group 2
|
9.651163
|
|
92
|
51
|
540
|
202
|
3
|
5.154e+11
|
Block Group 3
|
18.833333
|
|
108
|
51
|
540
|
202
|
1
|
5.154e+11
|
Block Group 1
|
8.282908
|
cville_blkgps$tree_can <- cvl_blkgps$tree_can
pal <- colorNumeric("plasma", reverse = TRUE, domain = cvl_blkgps$tree_can)
cvl_blkgps$tree_can <- as.numeric(cvl_blkgps$tree_can)
cville_blkgps$COUNTYFP <- as.numeric(cville_blkgps$COUNTYFP)
m <- leaflet()%>%
addTiles()%>%
addPolygons(data = cville_blkgps,
fillColor = ~pal(tree_can),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(weight = 2, fillOpacity = 0.8, bringToFront = T),
popup = paste0("FIPS Code: ", cvl_blkgps$GEOID, "<br>",
"Tree Canopy: ", cvl_blkgps$tree_can,"<br>",
"Blockgroup: ", cvl_blkgps$NAMELSAD)) %>%
addLegend("bottomright", pal = pal, values = cvl_blkgps$tree_can,
title = "Tree Canopy", opacity = 0.7)
m
Blocks
cvl_blocks <- read.csv("Cville_Tree/nlcd_tree_cville_blocks.csv")%>%
as.data.frame()
kable(head(cvl_blocks))
|
STATEFP10
|
COUNTYFP10
|
TRACTCE10
|
BLOCKCE10
|
GEOID10
|
NAME10
|
tree_can
|
|
51
|
3
|
10401
|
3025
|
5.100301e+14
|
Block 3025
|
74.666667
|
|
51
|
3
|
10401
|
3004
|
5.100301e+14
|
Block 3004
|
0.000000
|
|
51
|
3
|
10401
|
3005
|
5.100301e+14
|
Block 3005
|
25.907297
|
|
51
|
3
|
10401
|
3021
|
5.100301e+14
|
Block 3021
|
46.275294
|
|
51
|
3
|
10401
|
3016
|
5.100301e+14
|
Block 3016
|
44.680297
|
|
51
|
3
|
10401
|
3019
|
5.100301e+14
|
Block 3019
|
7.082353
|
# mean_block <- aggregate(cvl_blocks$tree_can, by = list(cvl_blocks$BLOCKCE10), FUN=mean)%>%
# round(3)
# colnames(mean_block) <- c("Blocks", "Tree_can")
cville_blocks$tree_can <- cvl_blocks$tree_can
pal <- colorNumeric("plasma", reverse = TRUE, domain = cvl_blocks$tree_can)
cvl_blocks$tree_can <- as.numeric(cvl_blocks$tree_can)
cville_blocks$COUNTYFP <- as.numeric(cville_blocks$COUNTYFP)
m <- leaflet()%>%
addTiles()%>%
addPolygons(data = cville_blocks,
fillColor = ~pal(tree_can),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(weight = 2, fillOpacity = 0.8, bringToFront = T),
popup = paste0("FIPS Code: ", cvl_blocks$GEOID, "<br>",
"Tree Canopy: ", cvl_blocks$tree_can,"<br>",
"Block: ", cvl_blocks$NAME10)) %>%
addLegend("bottomright", pal = pal, values = cvl_blocks$tree_can,
title = "Tree Canopy", opacity = 0.7)
m